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Quantum walks with long-range steps R ' {R being the distance between sites) on a discrete line behave in 
similar ways for all 7 > 2. This is in contrast to classical random walks, which for 7 > 3 belong to a different 
universality class than for 7 < 3. We show that the average probabilities to be at the initial site after time t as 
well as the mean square displacements are of the same functional form for quantum walks with 7 = 2, 4, and 
with nearest neighbor steps. We interpolate this result to arbitrary 7 > 2. 
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PACS numbers: 05.60.Gg, ()5.60.Cd, 03.67.-a, 71.35.-y 

One-dimensional models are not only prime toys for the- 
oretical physicists but also allow for deep physical insights. 
For instance, in solid state physics lattice models describe the 
behavior of metals quite acurately 1 1, 2]. Over the years these 
models have been refined and augmented to address different 
phenomena, such as the dynamics of atoms in optical lattices 
and the Anderson localization in systems with energetic dis- 
order [3]. Classical one-dimensional models allow to address 
various aspects of normal and anomalous diffusion |4]. 

The simplest model describing a particle moving on a reg- 
ular structure assumes only jumps from one position j to its 
nearest neighbors (NN) j ± 1. The tight-binding approxima- 
tion for such systems is equivalent to the so-called continuous- 
time quantum walks (CTQW), which model quantum dy- 
namics of excitations on networks jsl |6l I?!]. Recently, there 
have been several experimental proposals addressing CTQW 
in various types of systems, ranging from microwave cavities 
ySl], waveguide aiTays |9], atoms in optical lattices jld lllll . 
or structured clouds of Rydberg atoms [12]. A large class 
of these systems do not show NN steps. Consider, for in- 
stance, a chain of clouds of Rydberg atoms where each cloud 
can contain only one excitable atom due to the dipole block- 
ade lll2[ Il3n . The excited atoms of different clouds interact 
via long-range couplings decaying as R^^, where R is the 
distance between different clouds. 

The dynamics of classical excitations can be efficiently de- 
scribed by continuous-time random walks lll4il . Here, is has 
been shown that CTRW in one dimension with step lengths 
decaying as it!~^ belong only to the same universality class if 
7 > 3. Those CTRW show normal diffusion, whereas CTRW 
with 7 < 3 show anomalous diffusion as, e.g.. Levy flights. 
The reason is that the for 7 < 3 the second moment of the 
step-length distribution (i?^) diverges lll2lll5ll . 

In the following we will consider in one dimension the de- 
pendence of the dynamics of excitations on the range of the 
step length. We restrict ourselves to the extensive cases, i.e., 
we explicitly exclude ultra-long range interactions, where the 
exponent 7 of the decay of the step length is smaller than 
the dimension {-f < d, j — d is the marginal case); thus 
we take here 7 > 2. The effect of ultra-long range inter- 
actions on the thermodynamics and dynamics of regular one- 
dimensional lattices has been studied numerically before lll6ll . 

Our analysis is based on the density of states (DOS) of the 



coiTesponding Hamiltonian. The DOS contains the essential 
information about the system and allows to calculate various 
dynamical quantities, such as the probability to be at time t at 
the initially excited site. 

The coherent dynamics of excitons on a graph of connected 
nodes is modeled by CTQW, which follows by identifying the 
Hamiltonian H of the system with the CTRW transfer matrix 
T, i.e., H = -T; see e.g. [5, 6] (we wiU set /i ee 1 in the fol- 
lowing). For NN step lengths and identical transfer rates, T is 
related to the connectivity matrix A of the graph by T = —A. 
In the following, we will consider one-dimensional networks 
with periodic boundary conditions (i.e., a discrete ring). Here, 
when the interactions go as R^'^, with R = \k — j\ being the 
(on the ring minimal) distance between two nodes j and k, the 
Hamiltonian has the following structure: 



H, 



N flmax 



2. 



1 R=l 



R- 



2\n) {n\ -\n- R) {n\ - \n + R) {n\ 



(1) 
where i?niax is a cut-off for finite systems. Note, that in the 
infinite system limit we first take N ^ 00 before taking also 
^max -^ 00. For the cases considered here, namely 7 > 2 and 
N of the order of a few hundred nodes, a resonable cut-off is 
^max = N /2, which is also the largest distance between two 
nodes on the discrete ring. In this way, to each pair of sites a 
single (minimal) distance and a unique interaction is assigned. 
The states |j) associated with excitons localized at the 
nodes j {j = 1, . . . , N) form a complete, orthonormal ba- 
sis set (COBS) of the whole accessible Hilbert space, i.e., 
{k\j) = Skj and J^k 1^ )( ^1 = 1- ^^ general, the transition 
amplitudes from state \j) to state \k) during t and the corre- 



sponding probabilities read aZ (t) = 






(t) ^ 



^S 



iH.yi)|j) and 



{t)\ , respectively. 
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In the classical CTRW 



case the transition probabilities obey a master equation and 
can be expressed as p^^^(t) — (fc| exp(T-yi)|j) ||5l|6|]. 

For all 7, the time independent Schrodinger equation 
H-y|$e) = E^{9)\^g) is diagonalized by Bloch states |$e) = 
N^^^^ S7=i 6xp(i6'j)|j). One obtains the eigenvalues 



E^{9) = J2 ^"'^ [2 - 2 cos(6'i?)] . 



(2) 



R=l 



In the limit N ^ oo the 6'-values are quasi-continuous. 
Then, the density of states (DOS) p-y{E) is obtained by in- 
verting Eq. ^ and taking the derivative with respect to Ej. 
In the NN-case (7 = 00) only the first term in Eq. (|2]l, 
R = 1, contributes. From this we get the known DOS 
p^{E) = {ttV^E~~W)^\ For 7 = 2 we can ap- 
proximate the sum by letting iimax -^ cxd, which yields 
E2iO) = ne - 9^/2 (see Eq. 1.443.3 of fv^). By inverting 
this and assuming 9 to be continuous one obtains P2{E) ~ 
(7rV2\/7r2/2 — E) . In the intermediate range we have 
an analytic solution for 7 = 4, namely we have Ei{9) ~ 
6i4/24 - 776*3/6 + TT^e^Ve (see Eq. 1.443.6 of \m), which 

yields 01 pi{E) = [27r(2/3)i/S 



£;3/2 



'i;(7r2/V24) 

In order to interpolate between p2 [E) and p^o (E) to arbi- 
trary values of 7 e [2, 00] we assume the following general 
form for the DOS: 



p^{E) ^ 

with a e [0, 1] and /3 
the maximal energy, c^ 
a — and (3 = I (c2 
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[1,2]; c^ is a constant related to 

(^7,max)^^". Thus, for 7 = 2: 
77^/2). for 7 = 4: a = 1 and 



/3 = 3/2 (C4 = 7rV%/24), and for NN-walks: a 
P = 2 (coD = 4). Note that close to the band edge 9 
for small E, Eq. ^ can be approximated by 



PfiE) 



^-1/2 



1 (7 > 3) 



£;-"/2 a < 1 (2 < 7 < 3) 



1 and 
0, i.e.. 
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from which we observe the distinction between 7-values 
larger and smaller than three. For the band edge 9 ~ tt, 
i.e., E « -Bmaxj it is straightforward to show that p{E) ^ 
(^max - -E)-i/2 for all 7 > 2. The behavior of Eq. Q for 
small E is in line with previous studies lll9[l . in which the 
DOS goes as p{E) - E", where j/ ^ -1/2 for 7 > 3 and 
1/ — —(7 — 2)/(7 — 1) for 2 < 7 < 3. Starting from the two 
limiting cases 7 = 2 and 7 = 00 and supported by the 7 = 4 
case, Eq. ^ appears as a natural candidate for a generalized 
DOS. 

Figure[T]shows a comparison of the DOS obtained from the 
numerical diagonalization of H-,, for N = 10000 with 7 = 2, 
3, 4, and 00 (solid black curves) with the exact expressions 
for P2{E), Pi{E), and p^{E), see above, as well as a fit for 
PsiE). The values of a and (3, extracted from fits to the nu- 
marical DOS for various values of 7, are given in the inset of 
Fig.[T] Clearly, for 7 > 4 we have a = 1, while f3 e]l, 2[. 
For 7 = 2, the values of a and f3 drop to a = and /5 = 1, 
respectively. 

CTRW with step widths distributed according to R^ belong 
to the same universality class for 7 > 3, the mean square dis- 
placement (MSD) going as (i?^) ^ t, i.e., showing normal 
diffusion, see e.g. Ill4tl . For 7 < 3 the second moment of the 
distribution diverges, which leads to a MSD showing anoma- 
lous diffusion. 

Another way to see this is using the average probability to 
be at the initial site at time t, p~{t). Classically one has a 




FIG. 1: (Color online) Density of states for a discrete ring with 
A^ = 10000 nodes with hopping parameters 7 = 2, 3, 4, and 00 
(solid black curves with symbols), obtained by numerically diago- 
nalizing H^. The dashed green curves show the analytic expressions 
for p2{E), pa{E), and poo(-E) as well as the fit for p-j,{E) [Eq. ©l 
given in the text. The inset shows the exponents a and /3 obtained by 
fitting the DOS for different values of 7 to Eq. ^. 



simple expression forp (t) ll20ll2l|] 
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which depends only on the eigenvalues but not on the eigen- 
vectors. In the quantum case, the corresponding expression is 

7f^ {t) = j^ Ylf=i ^jj (*)■ For the discrete ring, we get 



n^{t) = |a7Wl' = I^E exp[-i£;^(0)i] 



(6) 



which also depends only on the eigenvalues. Note that for 
more complex networks the right-hand-side of Eq. (|6]l is only 
a lower bound to nj{t) 12211 . In the continuum limit, Eqs. (|5]i 
and ^ can be written as 



P^{t) = I dEp^{E) exp{-Et), 



TT^it) 



dE p^{E) e:cp{-iEt) 



(7) 



(8) 



Having the DOS at hand the integrals in Eqs. (|7]i and (O 
can be calculated - at least asymptotically - for large t. In the 
classical case Eq. (|7]i will be dominated by small values of E 
when t becomes large, see Eq. (|4]l. From the DOS we obtain 



Mt) 



t-V2 a = 1 
^a/2-1 Q, ^ 2_ 
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Quantum mechanically, some care is in order. Here, the 
assumption that Woo (t) will be dominated by small values of E 
for large t is not applicable, due to the oscillating exponential 
in Eq. ^. For the NN-case we know that 7foo(t) ~ t~^, see 
for instance [22]. Considering now the other limiting case 



2 we have 



^2{t) = 



7rV2 



dE e^p{-iEt)/TTV2^/Tr'^/2 - E 



TT^/2 



de exp(— iet) /v27r^ 



t-\ (10) 



where we substituted e = 7r^/2 — £'. Note that for i :» 1, 
Eq. ( [Tol l approaches 7f2(t) ~ (27rt)^^. Thus, the dependence 
of 7r2(i) on t is the same as for 7foo(i). This suggests that for 
all one-dimensional lattices with extensive (7 > 2) interac- 
tions the long time dynamics of the excitations is similar, no 
matter how long- or short-range the step lengths are. This is in 
contrast to the classical case, where only CTRW with 7 > 3 
belong to the same universality class. 
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FIG. 2: (Color online) (a) Classical p^(t) and (b) corresponding 
MSD; (c) quantum mechanical 7f-y(i) and (d) corresponding MSD 
(right) for a discrete ring with A*' = 10000 nodes with 7 = 2, 3, 4, 
and cx). 

To test this we calculated numerically for a discrete ring of 
N = 10000 nodes p^(i) and 7f^(i) for different 7; the results 
are shown in Fig.|2l Clearly, p^(t) changes when increasing 
the step width from NN steps (7 — co) to long-range steps 
distributed as R^"^ (7 = 2), see Fig. I^a). While f^{t) for 
7 > 3 decays as i^^/^, the power law changes to t^^ for 
7 = 2. In contrast, the decay of the maxima of the quantum 
return probability 7f^ (t) follows tr^ for all 7, Fig.|2c). Long- 
range steps lead only to a damping of the oscillations and to 
an earlier interference once the excitation has propagated half 
around the ring. 

The classical and quantum MSD corroborate these findings, 
see Figs.Eb) andEd). Now, the MSD for CTRW/CTQW on 
the discrete ring with initial site j are given by 



(i??W)c 



1 ^ 

-y 



MSD. For 2 < 7 < 3 the MSD starts to diverge, which in the 
case of finite networks is reflected in the fact that the MSD is 
of the order of N already for very short times. Increasing 7 to 
values larger than 3 leads to the expected diffusive behavior 
(i?^(t))ci ~ t for all 7 > 3. The quantum MSD, on the other 
hand, do not diverge for all 7-values considered here. All step 



/qm 



t^ 



lengths lead to the same qualitative behavior, (i?^ (t)) 

Figure |2] also shows that the MSD can be related to V-^ify 
and7f-y(i), through: 
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(12) 



for 7 > 3 in the classical and 7 > 2 in the quantal case. This 
generalizes previous (classical) results, obtained for regular 
networks with NN-steps 1201 . to the quantum case and to long- 
range steps. 

We can now underline our results by analytically evaluat- 
ing 7f-y(t}^[Eq. (|6]l] using the stationary phase approximation 
(SPA) 12311 . We expect in general Ej{9) to be a smooth real- 
valued function on the interval 6 E [0, 27r[. For large N, we 
write a^ {t) [Eq. (|6]l] in the integral form 

r2iT 

de exp{iE^{0)t). 



"^^'^ = ^i 



The SPA asserts now that the main contribution to this in- 
tegral comes from those points where E^{6) is stationary 
[dE^{0)/d0 = E'^{9) = 0]. If there is only one point 6*0 
for which E'^ido) = and d^E^i0)/d0^\g, = E!^{9o) ^ 
one gets, see [23J, 

a^{t) « , /_ exp (i{tE^{0o) + ^sgn[E'^{0o)]}) 



such that 



2TTt\E''{0o)\ 



,{t) = \a^{t)\' 



2Trt\E'^{ 



r 



(13) 



(14) 



For the infinite one-dimensional regular network [see Eq.©, 
where N ^ 00 and R^ax — > 00] and for 7 = 2, E2{0) (see 
above) has only one stationary point in e [0, 27r[, namely 
6*0 = TT. Then E2{tt) = tt^S and E'^{-k) = -1, leading to 
7^2(0 ~ (27ri)^^, which does not show any oscillations and 
coincides with the long time limit of Eq. ( fTol i. 

For J > 2, E^ (0) [see Eq. (|2j] has two stationary points 
in the interval £ [0, 27r[, namely 0o = and 0q = 
TT. Then a^{t) is approximately given by the sum of the 
contributions [each being of the form given in Eq. ( fT3] ll 
of the two stationary points. One easily verifies from 
E!^{0) = 2Y^j^cos{0R)/R''-^\hatsgn[E'^{O)] = land 
sgn[_E''(7r)] = — 1. Consequently, we obtain 



1 



ik-ji^n-^'it), (11) Tf^w = i«7Wr«^ 
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where 7'^;^^^(i) 



,(7), 



A^)i 



p^^]{t) for CTRW and ^^^(i) = n)7>(t) 
for CTQW. Now, decreasig 7 has huge effects on the classical 



+ 



2cos{t[E^{0) - E^{tt)] + 7r/2}\ 



^\E^{0)Ei;i7T)\ 



(15) 



The results for the infinite system and arbitrary 7 > 2 are 
readily obtained: for 6*0 = we have E^{0) = for all 7 
and E!^{0) = 2C(7 - 2), where C(7}_^ E^=i R^'^ is the 



Riemann zeta function, Eq. 23.2.1 of BIST . For 6q ~ tt we get 

E^in) = S^,„,axand^;'(7r) = 2r;(7-2) = (2-24-^)C(7- 
2), where 77(7) = Ek=i(-1)^"^-R"^' Eq. 23.2.19 of ifljl . 
Hence 



7r^(i) 



-I- 

47ri I |C(7 
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2)1 b7(7-2)| 
2cos[t^^(7r)+7r/2] 1 
VlC(7-2)r](7-llTj 



(16) 



For 7 = 3, this yields TTs (i) w [27rln(2)t]^^, which also does 
not show any oscillations. Comparing Eq. ( fTSI l for 7 = cx) 
to the long-time of the exact solution I.2Z1 . we have 7foo(i) ~ 
[2 - 2cos(4t + 7r/2)]/(27rt) = sin2(2i + 7r/4)/(7ri), which 
is exactly the asymptotic expansion of 7foo(i) = | Jo(2t)p w 
sin^(2i + Tr/ 4)/{nt) , where J,„(2t) is the Bessel function of 
the first kind iflslB. 
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FIG. 3: (Color online) Comparison of 7f-y(i) obtained from exact 
diagonalization (solid black curves) and the SPA (dashed red curves) 
for 7 = 2, 3, 4, and 00. 

Figure [3] shows comparisons of 7r-y(t) obtained from the 
exact diagonalization of H^ (solid black curves) to the SPA 
[Eq. ( fTSI ), dashed red curves]. Clearly, the oscillations de- 
crease with decreasing 7. For 7 = 2 and large t, the os- 
cillations of the exact 7f2(t) have practically vanished and 
7F2(t) « (27rt)~^. The SPA for 7 = 3 still shows oscillations 
because E''{Q) < 00. Increasing 7 further leads to an even 
better agreement of the SPA with the numerically evaluated 
decay. 

In conclusion, we have analyzed the quantum dynamics of 
excitations on discrete rings under long-range step lengths. 



distributed according to R~''. For specific cases, we calcu- 
lated the DOS analytically and interpolated to arbitrary step 
length ranges. The analytically obtained DOS enabled us to 
analytically calculate the average probability to be at the ini- 
tial site at t, which we related to the mean square displace- 
ment at time t. The classical MSD show that only CTRW 
with 7 > 3 belong to the same universality class, displaying 
normal diffusion. In contrast, the quantal MSD increase as P 
for all extensive cases, 7 > 2. Analytic calculations of the 
probability to be at the initial node within the stationary phase 
approximations confirm these findings. 
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